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Abstract 

Propagation of a Brillouin precursor in a Lorentz dispersive me- 
dium is considered. The precursor is excited by a sine modulated 
initial signal, with its envelope described by a hyperbolic tangent func- 
tion. The purpose of the paper is to show how the rate of growth of 
the initial signal affects the form of the Brillouin precursor. Uniform 
asymptotic approach, pertinent to coalescing saddle points, is applied 
in the analysis. The results are illustrated with numerical examples. 
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1 Introduction 

It is well known that if very fast, rapidly oscillating signals propagate in a 
real medium, they undergo the dispersion phenomenon. Various frequency 
components of a signal propagate with different phase velocities, and they are 
differently dumped. As a result, the shape of the signal is distorted during 
propagation the signal in the medium. Naturally, this phenomenon is prac- 
tically important only at very short times and very high frequencies (of the 
order of 10^^ Hz and above in the assumed model). In now classical works 
Sommerfeld [1] and Brillouin [2, 3] have shown that in the Lorentz model 
of a dispersive medium, apart of the main signal two small precursors are 
formed. In the asymptotic description of the total field these precursors are 
interpreted as contributions to the field resulting from two different pairs of 
saddle points. For the Sommerfeld precursor pertinent simple saddle points 
vary outside some disc in a complex frequency plane. As the space-time 
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coordinate ^, to be defined later, takes the initial value equal unity, those 
points merge at infinity to form one saddle point of infinite order. As grows 
up to infinity, they separate into two simple saddle points that move sym- 
metrically with respect to the imaginary axis towards corresponding branch 
points located in the left and the right half-plane, respectively. In the case of 
Brillouin precursor, two other simple saddle points vary inside a smaller disc. 
As the coordinate 9 grows from unity, they move toward each other along 
the imaginary axis, coalesce into one saddle point of the second order on the 
axis, and then again split into simple saddle points that depart from the axis 
and move, symmetrically with respect to this axis, towards corresponding 
branch points in the left and the right half-plane, respectively. The location 
of the saddle points affects local oscillations and dumping of the precursor. It 
depends on the space-time coordinate 6 and is governed by the saddle point 
equation. 

In this paper we confine our attention to the Brillouin precursor, also 
called a second precursor (as opposed to the first, Sommerfeld precursor). 
Fundamental work on this precursor is due to Brillouin [2, 3]. Because of limi- 
tations of asymptotic methods then available (now referred to as non-uniform 
methods), Brillouin could not correctly describe the precursor's dynamics for 
values of 6 corresponding to the coalescence of simple saddle points into one 
saddle point of a higher order. With the development of advanced, uniform 
asymptotic techniques, complete description of the precursor now got feasi- 
ble (Kelbert and Sazonov [4], and Oughstun and Sherman [5]). In the latter 
monograph, in addition to the delta function pulse, the unit step-function 
modulated signal and the rectangular modulated signal, the authors also 
studied an initial signal with finite rate of growth. In their model, however, 
the envelope of the initial signal is described by everywhere smooth function 
of time, tending to zero as time goes to minus infinity ([5], Sec. 4.3.4). In 
the present paper we consider more realistic excitation which is caused by 
an abruptly switched modulated sine signal, vanishing identically for time 
i < and being non-zero for i > 0. At i = the derivative of the sig- 
nal's envelope suffers a step discontinuity. As t increases, the envelope grows 
with a finite speed, asymptotically tending to its maximum value. In the 
following sections we construct uniform asymptotic representation for the 
Brillouin precursor resulting from this sort of excitation, and show how the 
speed of growth in the initial signal affects the form of the precursor. We 
also illustrate the results with numerical examples. 



2 



2 Formulation of the problem 



We consider a one dimensional electromagnetic problem of propagation in a 
Lorentz medium. The medium is characterized by the frequency-dependent 
complex index of refraction 

/ 62 \ 1/2 

n(u;) = 1 - — ^ — , (1) 

where b is so called plasma frequency of the medium, 5 is a damping constant 

and ujq is a characteristic frequency. 

Any electromagnetic field in the medium satisfies the Maxwell equations 

_ , ldH(r,t) ^ _ , ldE(r,t) 

D(v, t) = f £(i - T)E(t, T)dT, B(v, t) = nH{v, t), 

where e{t) is a real function and /i. is a real constant (hereafter assumed to 
be equal 1). By Fourier transforming the equations with respect to t and 
assuming that the fields depend on one spatial coordinate z only, we obtain 
the following equations for transforms of the respected fields 

z X Wz^uj) = —8[z,uj), z X t[z,uj) = H{Z,U!), 

c c 

where z is the unit vector directed along z-axis and e(a;) = n'^{u!)/{c^ii) is 
the Fourier transform of e{t). It then follows that z, S and Ti are mutually 
perpendicular. Moreover, if S is known then Ti. is also known, and vice versa. 
It is also true for the electromagnetic field components, which are the inverse 
Fourier transforms of S and 7i. Therefore, the knowledge of the electric 
(magnetic) field is sufficient to determine the full electromagnetic field. To 
make the calculations as simple as possible, it is advisable that the x (or y) 
axis be directed to coincide with the electric or magnetic field. 

Assume that in the plane 2; = an electromagnetic signal is turned on at 
the moment t = 0. For t > it oscillates with a fixed frequency lJc and its 
envelope is described by a hyperbolic tangent function. Suppose the selected 
Cartesian component (say x-component) of one of these fields in the plane 
2; = is given by 

. . . _ f t <0 , . 

~ \ tanhptsincu^t t > 0, ^ ^ 

The parameter P determines how fast the envelope of the signal grows. 
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This initial electromagnetic disturbance excites a signal A{z, t) outside the 
plane z = 0. In what follows we will be interested in the field propagating in 
the half-space z > 0. The problem under investigation can be classified as a 
mixed, initial-boundary value problem for the Maxwell equations. 

The exact solution for this specific form of the initial signal A{0,t) is 
described by the contour integral [6] 



J c 

defined in the complex frequency plane uj. Here, 



(3) 



1 

47r 




2(5 , 



+ 



UJ — UJr 




i{uj -\- ujc) 



UJ -\- UJc 



the complex phase function 0(ti;, 6) is given by 



C r7 



[k{ij)z — ujt] — iuj[n{uj) — 6], 



and B{s) is the beta function [7] defined via the psi function as 



Bis) 



s + l 



(4) 



(5) 



(6) 



B{s) is a Fourier transform of the initial signal envelope tanh I3t. The dimen- 
sionless parameter 

p) 

defines a space-time point (z, t) in the field, and c is the speed of light in 
vacuum. The contour C is the line uj — uj' + ia, where a is a constant greater 
than the abscissa of absolute convergence for the function in square brackets 
in (4) and uj' ranges from negative to positive infinity. 

Our goal is twofold. First, we shall seek an asymptotic formula for the 
second (Brillouin) precursor that results from the excitation A(0, t). In other 
words, we shall find near saddle points contribution to the uniform asymp- 
totic expansion of the total field A{z,t). Second, we shall examine how the 
speed parameter /3 in (2) affects the form of the Brillouin precursor. 



3 Asymptotic representation for the second 
precursor 

Our derivation of the asymptotic formula for the Brillouin precursor is based 
on the technique developed by Chester et al. [8] for two simple saddle points 
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coalescing into one saddle point of the second order. The technique is also 
conveniently described in [9] and [10]. 

The locations in the complex cj-plane of the saddle points in (3) are 
determined from the saddle point equation 

n{u)+un'{u) -e ^0. (8) 

At these points the first derivative (j)^^{uj,6) of the phase function vanishes. 
We are interested in the near saddle points, varying in the domain < 



ujI + 5'^. As 9 increases from 1 to a value denoted by 9i, the near saddle 
points cui and CU2 approach each other along the imaginary axis from below 
and from above, respectively ([5]). They coalesce to form a second order 
saddle point at = ^1. Finally, as 6 tends to infinity they depart from the 
axis and symmetrically approach the points cui 2 = ±^0 ~ in the right 
and in the left complex cu half plane, respectively. If n{uj) is eliminated from 
(8) then the equation can be represented in the form of an eighth degree 
polynomial in u on its left hand side, and zero on its right hand side. It does 
not seem to be possible to solve the equation exactly. In what follows we shall 
employ the solution to (8) which was obtained numerically. Alternatively, a 
simple approximate solution found in [11] could be used here at the expense 
of accuracy in resulting numerical examples. 

The first step in the procedure is to change the integration variable t in 
(3) to a new variable s, so that the map s{uj) in some disk D containing the 
saddle points uj± (but not any other saddle points) is conformal and at the 
same time the exponent takes the simplest, polynomial form 

0(u;,^) = p + 7's-y = t(s,^). (9) 

Notice that r(s, 6) has two simple saddle points s = ±7 that can coalesce 
into one saddle point s = of the second order, corresponding to uj = LUg. 
Prom 

we infer that for s{uj) to be conformal, 5 = 7 should correspond to a; = 
and s — —7 should correspond to a; = a;2. Then, 



^^±^) = \/;a^' ^^^^ 

V 4> ('^1,2) 

where (j)" {^01^2) is a short notation for (f)" (uji^2:9). In case the saddle points 
cu — UJ12 merge to form one saddle point of the second order cu — cUs, one has 
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0" (ujs) — 0, and the relevant formula for (jg is 



a;(0) 



nl/3 



;i2) 



By using correspondence 001^2 ^ =1=7 in (9) one finds that 7^ and p are equal 
to 

^ = (13) 

+ (14) 

The equation (13) for 7 has three complex roots. Only one root corresponds 

to a regular branch of the transformation (9) leading to the conformal map 
s{uj). To find the proper value of 7 we first note from (12) that argcj(O) 
can take one of the three values: tt/G, Stt/G or — 7r/2 corresponding to three 
different branches of the transformation (9). It can be readily verified that 
for 9 < 9s both 0(a;i) and 0(a;2) are real valued and 0(a;i) > 0, 0(a;2) < 
0. Then it follows from (13) that 7^ > 0. On the other hand, if ^ > 
then uji = —002-^ This implies that in the present case 0(c<Ji) = [0(c(;2)]*, 
and similarly (cui) = [0 {uj2)]*- It is now seen that RHS of (9) equals 
—2i Im (t){uJ2), where Im 0(c<;2) < 0. Hence for 9 > 9s, arg7^ = 7r/2. We 
now take advantage of the fact that a;(s) as given by (11) tends in the limit 
to (12) as s — > 0. Because arg0(c<Ji) = and arg0(cj2) = tt for 9 < 9s, 
and arg(/)(cL;i) — —n/2 and arg(6(cj2) — >■ vr/2 as 6* ^ 9f, we conclude that 
arg7 = 6^ < and -7r/2 for 9 > 9s, i.e. 

7 = [||0(^i)- 0(^2)1]^/^ (15) 

where a = 0, — 7r/2 if ^ < or 9 > 9s, respectively. 

With the new variable of integration the integral (3) can be written down 
in the form 

A(z,t)= [ Gis,9)e^^^''^Us + 8, (16) 

where X — z/c and 

G{s,9)=g[u:{s)]^{s). (17) 

The contour Ci is an infinite arc in the left s complex half-plane, symmetrical 
with respect to the real axis, running upwards and having rays determined by 
the angles — i27r/3 and i27r/3 as its asymptotes. The domain D is the image 
of D under (9). The term £, standing for the integral of G{s, 9) exp At(s, 9) 

^The star denotes complex conjugate. 
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defined over the parts of Ci outside D, is exponentially smaller than A{z, t) 
itself. 

We now represent 9) in the canonical form 



G{s, e)^Co + CiS + {s^ - -f^)H{s, 0). 



(18) 



Provided the function H{s, 6) is regular, the last term in (18) vanishes at 
the saddle points s — ±7, and its contribution to the asymptotic expansion 
is smaller than that from the first two terms. Indeed, it can be shown that 
integration by parts of the last term leads to an integral of similar form as 
(11) multiplied by A^^. 

To determine cq and ci we substitute s = ±7 in (18) and thus find 



Co 



Cl 



G{j,9)-G{-j,9) 
27 



(19) 



(20) 



By using (18) and (9) in (11), and extending the integration contour in the 
resulting integrals to Ci, we find that the leading term of the asymptotic 
expansion of A{z, t) as A ^ 00 is given by 

A{z,t) - 2me^P^'^ (X-'^^co{9)Ai[X^/^'y{9f] + X-^/^ci{9)Ai'[X^/^'y{9f]) . 

(21) 

It is defined through the Airy function and its derivative, as given by ([10]) 



Ai(x) = —f e'^-'''^ds Ai'(x)^—f se'^^-^'/^ds. 
^ ' 2mJci 2mJci 



(22) 



Plots of both functions for real x are shown in Fig. 1 . 

The expansion holds for any 7(6'), including 7 = 0. This special case 
corresponds to coalescing of the two simple saddle points r = ±7 into one 
saddle point of the second order. In other words the expansion is uniform 
in 7, and hence in 9. It is seen that for 7 = 0, i.e. for 9 = 9s, the algebraic 
order of A(z,t) in A is A~^/^. This behavior is characteristic of an integral 
with a saddle point of the second order. 

For 7 well separated from zero the Airy function and its derivative can 
be replaced by their asymptotic expansions ([12]) 



Ai(a;) 



.3/2 



~ e 



Ai'{x) 



x3/2 



~ e 



2V7rxV4 

x^ 



+ 0(x-3/4) 

+ o(x-^/^) 



(23) 
(24) 
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Figure 1: Plots of Ai and Ai' against a real argument x. 



as X — > oo, and 



Ai{x) 



A(-x)V4 



sm 



(25) 



+ 



-X) 



-3/2 



]}. 



Ai'(a;) ~ 
+ 



-X 



)l/4 



-x)-'/']} 



COS 



{l + 0[{-x)-']) 



(26) 



as x —oo. By using these expansions in (21) we arrive at the following 
non-uniform asymptotic representation of the precursor 



1/2 



A{z,t) ~ e 



(27) 



if ^ < 9s, and 



' -27r \ 



1/2 



5f(u;i) +e 



A0(a)2) 



^ -27r \ 



1/2 



^(^2) (28) 



if ^ > e^. 

We see from the above formulas that for 6 sufficiently distant from 9s (for 
Brillouin's choice of medium parameters 9s ~ 1.5027), the representation 
(21) reduces to a simple saddle point contribution from a;2 if ^ < 9s, and to 
a sum of simple saddle point contributions from cui and CU2 H 9 > 9s- In this 
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manner it is confirmed that the saddle point uji does not contribute when 

9 < Off. This is a direct consequence of the fact that the original contour 
of integration in (3) cannot be deformed to a descent path from imaginary 
cu = cui. The algebraic order of A{z.t) in A is now A^^/^ because in this 
case separate simple saddle points contribute to the expansion. From (27) 
and (28) it is also seen that these formulas are non-apphcable at ^ = (i.e. 
7 = 0), where (1^1,2) = 0. On the other hand the uniform expansion (21) 
remains valid for any 6 (and 7). In particular it provides a smooth transition 
between the cases of small and large |7|. 

If 6* > 9s, then it can be readily seen that g{uJi) = g*{uj2)i (pi^i) = 4'*{^2)-i 
and similarly (a;i) = ('^2)- In this case (28) can be written down in a 
more compact form 



A{z, t) ~ 2Re 



-2ti 



1/2 



\(t)"{uJ2), 



(29) 



Fig. 2 shows the dynamics of the Brillouin precursor in a Lorentz medium 
as given by its uniform representation (21) and non-uniform ones (27) and 
(28). Throughout this work the Brillouin's choice of medium parameters 



h = X cuo = 4.0 X 10^^s-\ 5 = 0.28x10^^^8"^ (30) 

is assumed. 



A(z,t) 



0.15 



0.1 



0.05 



-0.05 




Figure 2: Uniform (solid line) and non-uniform (dashed line) representa- 
tion of the Brillouin precursor in a Lorentz medium described by Brillouin's 
parameters. Here, (3 = l.Ox lO^^s'S = 2.5x lO^^s"^ and A = 3.0x lO'^^s. 



For 9 < 9s the function 7(6*) takes positive values and the precursor is 
described by a monotonically changing function. Adversely, for 9 > 9s the 
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argument in both functions takes negative values which leads to oscillatory 
behavior of the precursor. This reflects the behavior of both Airy function 
and its derivative for positive and negative values of their argument (see 
Fig 1). 

4 Dependance of the precursor on the rate 
parameter (3 

A(z,t) 
0.0006 
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. 0002 



-0.0002 
-0.0004 
-0.0006 



A(z,t) 




e 



Figure 3: Dynamic behavior of the Brillouin precursor in the Lorentz medium 
described by Brillouin' s parameters obtained for: (top) j3 = 2.0 x lO^'^s"^, and 
(bottom) /3 = 2.0 X lO^^s-^ Here, cUc = 2.5 x lO^S"^ and A = 3.0 x lO-^^s. 

An important question arises on how the rate parameter /3 affects the 
form of the Brillouin precursor. 
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As the parameter /3 in (21) increases starting from relatively small values, 
the shape of the precursor remains virtually unchanged while its magnitude 
grows. This tendency is no longer valid if f3 enters a transitory interval. In 
that interval the shape of the precTirsor changes and its magnitude rapidly 
increases. Above transitory interval, further increase of (3 leaves the shape 
and the magnitude of the precursor virtually constant. The form of Brillouin 
precursor for (3 below (2.0 x lO^^s"^) and above (2.0 x lO^^s"^) the transitory 
interval is shown in Fig. 3. 

Explanation of this behavior lies in the properties of the coefficients Cq 
and Ci in (21), which are /3-dependent. The coefficients, in turn, determine 
the weight with which Airy function and its derivative contribute to the 
precursor. 




Figure 4: Plots of 2Txi\ ^f^CQ{9]j3) and 2-KiX '^^^Ci{9; (3) against the speed 
parameter (3 at 6 = 1.502. Here, Uc = 2.5 x lO^^s^^ and A = 3.0 x lO^^^s. 

First, consider the case of ^ < In Fig. 4 the coefficients of, respectively, 
Ai and Ai' in the parentheses in (21) multiplied by 27ri are plotted against 
(3. The value of 6 is chosen to be slightly below Og. For relatively small f3 the 
term proportional to Ai'(-) dominates over the term proportional to Ai(-). In 
this case the ratio of both terms remains unchanged in a wide interval of (3 
variation. The magnitude of the precursor increases with f3 growth up to the 
moment where the contribution from Ai(-) changes sign and rapidly grows 
until finally it settles down at a virtually constant level. At the same time the 
contribution from Ai'(-) decreases to another constant level and is very small 
compared to the other term. At this stage the shape and the magnitude 
of the precursor are approximately determined by the special form of the 
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function 

g{u;) = -L L\ (31) 

which appears in cq and ci, which is a hmiting case of g{u!) as /3 — > cxo, i.e. 
for the initial signal with a unit step function envelope. 

Now consider the case of ^ > ^5 where the precursor becomes oscillatory. 
The envelope of the oscillations can be conveniently approximated with the 
help of (29) by 

A{z, t; (5) ^ 2e^^^^^^'^^ 

provided |7(^)| is sufficiently large. 



-27r 



1/2 



^A0"(u;2), 



(32) 
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Figure 5: Dependence of the magnitude of the Brillouin precursor's envelope 
on (5. Calculated at different values of 6. Here, ujc — 2.5 x lO^^s"-*^ and 
A = 3.0 X 10-1^ 

In Fig. 5 the magnitude of the precursor envelope is plotted against the 
parameter /3 for different values of ^. It is seen again that after fast growth 
of the envelope magnitude at relatively small values of /?, which occurs with 
approximately the same rate for all 9, the magnitude reaches a saturation 
level at higher values of /3. Since the saturation appears earlier at larger 
values of 9, the precursor envelope has a tendency to become narrower with 
growing (5. Additionally, one observes that with growing (5 the first extremum 
moves towards larger values of 9. This is a direct consequence of the fact that 
the first extremum of the Airy function is shifted towards negative values of 
its argument as compared to the first extremum of the derivative of the Airy 
function. It has an additional effect on narrowing the precursor shape. 
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5 Conclusions 



In this paper we have derived the uniform and non-uniform asymptotic rep- 
resentations for the Brillouin precursor in a Lorentz medium, excited by an 
incident signal of finite rise time, and well defined, startup time. With the 
use of these representations we analyzed the effect of the speed parameter P 
on the form and magnitude of the precursor. The results obtained can be 
helpful e.g. in applications involving triggering devices that work with signal 
amplitudes close to the noise level. In this paper we did not consider the 
problem of smooth transition from Brillouin precursor to the main signal. 
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